Kinetic Properties Study of H Atom Abstraction by CH3Ȯ2 Radicals from Fuel Molecules with Different Functional Groups.

The detailed kinetic properties of hydrogen atom abstraction by methylperoxy (CH3Ȯ2) radicals from alkanes, alkenes, dienes, alkynes, ethers, and ketones are systematically studied in this work. Geometry optimization, frequency analysis, and zero-point energy corrections were performed for all species at the M06-2X/6-311++G(d,p) level of theory. The intrinsic reaction coordinate calculation was consistently performed to ensure that the transition state connects the correct reactants and products, and one-dimensional hindered rotor scanning results were performed at the M06-2X/6-31G level of theory. The single-point energies of all reactants, transition states, and products were obtained at the QCISD(T)/CBS level of theory. High-pressure-limit rate constants of 61 reaction channels were calculated using conventional transition state theory with asymmetric Eckart tunneling corrections over the temperature range of 298.15-2000 K. Reaction rate rules for H atom abstraction by CH3Ȯ2 radicals from fuel molecules with different functional groups are constructed, which can be used in the development of combustion models of these fuels and fuel types. In addition, the influence of the functional groups on the internal rotation of the hindered rotor is also discussed.

Experimental and Updated Kinetic Modeling Study of Neopentane Low Temperature Oxidation.

Neopentane is an ideal fuel model to study low-temperature oxidation chemistry. The significant discrepancies between experimental data and simulations using the existing neopentane models indicate that an updated study of neopentane oxidation is needed. In this work, neopentane oxidation experiments are carried out using two jet-stirred reactors (JSRs) at 1 atm, at a residence time of 3 s, and at three different equivalence ratios of 0.5, 0.9, and 1.62. Two different analytical methods (synchrotron vacuum ultraviolet photoionization mass spectrometry and gas chromatography) were used to investigate the species distributions. Numerous oxidation intermediates were detected and quantified, including acetone, 3,3-dimethyloxetane, methacrolein, isobutene, 2-methylpropanal, isobutyric acid, and peroxides, which are valuable for validating the kinetic model describing neopentane oxidation. In the model development, the pressure dependencies of the rate constants for the reaction classes Q̇OOH + O2 and Q̇OOH decompositions are considered. This addition improves the prediction of the low-temperature oxidation reactivity of neopentane. Another focus of model development is to improve the prediction of carboxylic acids formed during the low-temperature oxidation of neopentane. The detection and identification of isobutyric acid indicates the existence of the Korcek mechanism during neopentane oxidation. Regarding the formation of acetic acid, the reaction channels are considered to be initiated from the reactions of ȮH radical addition to acetaldehyde/acetone. This updated kinetic model is validated extensively against the experimental data in this work and various experimental data available in the literature, including ignition delay times (IDTs) from both shock tubes (STs) and rapid compression machines (RCMs) and JSR speciation data at high temperatures.

Pyridinium-Inspired Organocatalysts for Carbon Dioxide Fixation: A Density Functional Theory Inspection.

The current project aims to apply the virtues of minimalism to examine the catalytic ability of commercially organic compounds of small chemical structures to catalyze the coupling reaction between carbon dioxide and propylene oxide (PO) under mild conditions. The proposed catalysts are pyridinium iodide (A), 2-hydroxypyridinium iodide (B), and piperidinium iodide (C), where their structure is based on cooperative acidic and nucleophilic motifs. The quantum chemistry model, M062X-D3/def2-TZVP//M062X-D3/def2-SVPP, was used to understand the reaction mechanism and the catalytic performance. Since the coupling reaction was performed under excess PO, we proposed that PO serves as a reactant and solvent. Therefore, calculations were performed in gas and liquid phases for comparison. The findings indicated that the rate-determining step depends on the chemical structure of the catalyst and whether the phase is a gas or liquid phase. In general, modeling in the liquid phase produces potential energy surfaces of lower energy barriers. The noncovalent interactions reflect the role of hydrogen bonding in controlling the kinetic behavior of the coupling reaction. Based on the finding, catalyst A is the best candidate for transforming CO2 into cyclic carbonates.

High-Accuracy Heats of Formation for Alkane Oxidation: From Small to Large via the Automated CBH-ANL Method.

It is generally challenging to obtain high-accuracy predictions for the heat of formation for species with more than a handful of heavy atoms, such as those of importance in standard combustion mechanisms. To this end, we construct the CBH-ANL approach and illustrate that, for a set of 194 alkane oxidation species, it can be used to produce ΔHf(0 K) values with 2σ uncertainties of 0.2-0.5 kcal mol-1. This set includes the alkanes, hydroperoxides, and alkyl, peroxy, and hydroperoxyalkyl radicals for 17 representative hydrocarbon fuels containing up to 10 heavy atoms with various degrees of branching in the alkane backbone. The CBH-ANL approach, automated in the QTC and AutoMech software suites, builds balanced chemical equations for the calculation of ΔHf(0 K), in which the reference species may be up to five heavy atoms. The high-level ANL0 and ANL1 reference ΔHf(0 K) values are further refined for even the largest of these reference species with a novel laddering approach. We perform a comprehensive quantification of the uncertainties for both the individual reference species (the largest of which is 0.15 kcal mol-1) and the propagation of those uncertainties when used in the calculation of ΔHf(0 K) for the 194 target species. We examine the sensitivity of the predicted ΔHf(0 K) values to (i) electronic energies from various methods, including ωB97X-D/cc-pVTZ, B2PLYP-D3/cc-pVTZ, CCSD(T)-F12b/cc-pVDZ-F12//B2PLYP-D3/cc-pVTZ, and CCSD(T)-F12b/cc-pVTZ-F12//B2PLYP-D3/cc-pVTZ; (ii) the zero-point vibrational energies (ZPVEs), where we consider harmonic ZPVEs as well as two scaling-based estimates of the anharmonic ZPVEs, all implemented for both ωB97X-D/cc-pVTZ and B2PLYP-D3/cc-pVTZ calculations; (iii) the particular CBH-ANL scheme employed; and (iv) the procedure for choosing the reference conformer for the analyses. The discussion concludes with a summary of the estimated overall uncertainty in the predictions and a validation of the predictions for the alkane subset.

Hierarchical Study of the Reactions of Hydrogen Atoms with Alkenes: A Theoretical Study of the Reactions of Hydrogen Atoms with C2-C4 Alkenes.

The present study complements our previous studies of the reactions of hydrogen atoms with C5 alkene species including 1- and 2-pentene and the branched isomers (2-methyl-1-butene, 2-methyl-2-butene, and 3-methyl-1-butene), by studying the reactions of hydrogen atoms with C2-C4 alkenes (ethylene, propene, 1- and 2-butene, and isobutene). The aim of the current work is to develop a hierarchical set of rate constants for H atom addition reactions to C2-C5 alkenes, both linear and branched, which can be used in the development of chemical kinetic models. High-pressure limiting and pressure-dependent rate constants are calculated using the Rice-Ramsperger-Kassel-Marcus (RRKM) theory and a one-dimensional master equation (ME). Rate constant recommendations for H atom addition and abstraction reactions in addition to alkyl radical decomposition reactions are also proposed and provide a useful tool for use in mechanisms of larger alkenes for which calculations do not exist. Additionally, validation of our theoretical results with single-pulse shock-tube pyrolysis experiments is carried out. An improvement in species mole fraction predictions for alkene pyrolysis is observed, showing the relevance of the present study.

Theoretical Study of the Reaction of Hydrogen Atoms with Three Pentene Isomers: 2-Methyl-1-butene, 2-Methyl-2-butene, and 3-Methyl-1-butene.

This paper presents a comprehensive potential energy surface (PES) for hydrogen atom addition to and abstraction from 2-methyl-1-butene, 2-methyl-2-butene, and 3-methyl-1-butene and the subsequent ß-scission and H atom transfer reactions. Thermochemical parameters for species on the C5H11 potential energy surface (PES) were calculated as a function of temperature (298-2000 K), using a series of isodesmic reactions to determine the formation enthalpies. High-pressure limiting and pressure-dependent rate constants were calculated using Rice-Ramsperger-Kassel-Marcus theory with a one-dimensional master equation. A number of studies have highlighted the fact that C5 intermediate species play a role in polyaromatic hydrocarbon formation and that a fuel's chemical structure can be key in understanding the intermediate species formed during fuel decomposition. Rate constant recommendations for both H atom addition to, and H-atom abstraction by H atoms from, linear and branched alkenes have subsequently been proposed by incorporating our earlier work on 1- and 2-pentene, and these can be used in mechanisms of larger alkenes for which calculations do not exist. The current set of rate constants for the reactions of H atoms with both linear and branched C5 alkenes, including their chemically activated pathways, are the first available in the literature of any reasonable fidelity for combustion modeling and are important for gasoline mechanisms. Validation of our theoretical results with pyrolysis experiments of 2-methyl-1-butene, 2-methyl-2-butene, and 3-methyl-1-butene at 2 bar in a single pulse shock tube (SPST) were carried out, with satisfactory agreement observed.

Hindered rotor benchmarks for the transition states of free radical additions to unsaturated hydrocarbons.

The hindered internal rotors of 32 transition states (TSs) formed through four free radicals, namely methyl, vinyl, ethyl, methoxy (CH3, C2H3, C2H5, CH3) additions to acetylene, ethylene, allene, propyne, and propene (C2H2/C2H4/C3H4-a/C3H4-p/C3H6) are studied. To validate the uncertainties of rate constants that stem from the use of different electronic structure methods to treat hindered rotors, the rotations of the newly formed C-C and/or C-O rotors in the transition states are calculated using commonly used DFT methods (B3LYP, M06-2X, ωB97X-D and B2PLYP-D3 with two Pople basis sets (6-31+G(d,p), 6-311++G(d,p)) and cc-pVTZ). The hindrance potential energies V(χ) calculated using the M06-2X/6-311++G(d,p) method are benchmarked at the CCSD(T), CCSD(T)-F12, DLPNO-CCSD(T) levels of theory with cc-pVTZ-F12 and cc-pVXZ (X = T, Q) basis sets and are extrapolated to the complete basis set (CBS) limit. The DLPNO-CCSD(T)/CBS method is proven to reproduce the CCSD(T)/CBS energies within 0.5 kJ mol-1 and this method is selected as the benchmark for all of the rotors in this study. Rotational constants B(χ) are computed for each method based on the optimized geometries for the hindrance potential via the I(2,3) approximation. Thereafter, the V(χ) and B(χ) values are used to compute hindered internal rotation partition functions, QHR, as a function of temperature. The uncertainties in the V(χ), B(χ) and QHR calculations stem from the use of different DFT methods for the internal rotor treatment are discussed for these newly formed rotors. For rotors formed by + C2 alkenes/alkynes, the V(χ) and QHR values calculated using DFT methods are compared with the DLPNO-CCSD(T)/CBS results and analysed according to reaction types. Based on comparisons of the DFT methods with the benchmarking method, reliable DFT methods are recommended for the treatment of internal rotors for different reaction types considering both accuracy and computational cost. This work, to the authors' knowledge, is the first to systematically benchmark hindrance potentials which can be used to estimate uncertainties in theoretically derived rate constants arising from the choice of different electronic structure methods.

An ab Initio/Transition State Theory Study of the Reactions of C5H9 Species of Relevance to 1,3-Pentadiene, Part II: Pressure Dependent Rate Constants and Implications for Combustion Modeling.

The temperature- and pressure-dependence of rate constants for several radicals and unsaturated hydrocarbons reactions (1,3-C5H8/1,4-C5H8/cyC5H8 + H, C2H4 + C3H5-a, C3H6 + C2H3) are analyzed in this paper. The abstraction reactions of these systems are also calculated and compared with available literature data. C5H9 radicals can be produced via H atom addition reactions to the pentadiene isomers and cyclopentene, and also by H-atom abstraction reactions from 1- and 2-pentene and cyclopentane. Comprehensive C5H9 potential energy surface (PES) analyses and high-pressure limiting rate constants for related reactions have been explored in part I of this work ( J. Phys. Chem. A 2019, 123 (22), 9019-9052). In this work, a chemical kinetic model is constructed based on the computed thermochemistry and high-pressure limiting rate constants from part I, to further understand the chemistry of different C5H8 molecules. The most important channels for these addition reactions are discussed in the present work based on reaction pathway analyses. The dominant reaction pathways for these five systems are combined together to generate a simplified C5H9 PES including nine reactants, 25 transition states (TSs), and nine products. Spin-restricted single point energies are calculated for radicals and TSs on the simplified PES at the ROCCSD(T)/aug-cc-pVTZ level of theory with basis set corrections from MP2/aug-cc-pVXZ (where X = T and Q). Temperature- and pressure-dependent rate constants are calculated using RRKM theory with a Master Equation analysis, with restricted energies used for minima on the simplified C5H9 PES and unrestricted energies for other species, over a temperature range of 300-2000 K and in the pressure range 0.01-100 atm. The rate constants calculated are in good agreement with those in the literature. The chemical kinetic model is updated with pressure-dependent rate constants and is used to simulate the species concentration profiles for H atom addition to cyclopentane and cyclopentene. Through detailed analyses and comparisons, this model can reproduce the experimental measurements of species qualitatively and quantitatively with reasonably good agreement.

Ab Initio/Transition-State Theory Study of the Reactions of C5H9 Species of Relevance to 1,3-Pentadiene, Part I: Potential Energy Surfaces, Thermochemistry, and High-Pressure Limiting Rate Constants.

In this study, the reactions of C5H9 radicals are theoretically investigated, with a particular emphasis on hydrogen atom addition reactions to 1,3-pentadiene (C5H8) to form C5H9 radicals, although the subsequent isomerization and decomposition reactions of the C5H9 radicals are also of direct relevance to the radicals formed from the pyrolysis and oxidation of species including pentene and cyclopentane. Moreover, H-atom abstraction reactions by hydrogen atoms from 1,3-pentadiene are also investigated. The geometries and frequencies of 63 potential energy surface (PES) minima and 88 transition states are optimized at the ωB97XD/aug-cc-pVTZ level of theory. Spin-unrestricted open-shell single-point energies for all the species are calculated at the CCSD(T)/aug-cc-pVTZ level of theory with basis set corrections from MP2/aug-cc-pVXZ (where X = T and Q). A one-dimensional hindered rotor treatment is employed for torsional modes, with the M06-2X/6-311++G(D,P) method used to compute the potential energy as a function of the dihedral angle. The high-pressure limiting rate constants and the thermochemical properties for C5 species are calculated using the Master Equation System Solver (MESS) with conventional transition-state theory and comparisons made with existing available literature data. A hydrogen atom can add to the terminal carbon atom of 1,3-pentadiene to form the 2,4-C5H9 radical and/or the internal carbon atoms to form 2,5-C5H9, 1,4-C5H9, and 1,3-C5H9 radicals. Among the four entrance channels for H atom addition reactions, the formation of 2,4-C5H9 and 1,3-C5H9 radicals is more exothermic in comparison to the other C5H9 isomers (2,5-C5H9, 1,4-C5H9) because of the resonantly stabilized allylic structure. Consequently, the formation of the former is generally dominant in terms of barrier heights. H atom addition reactions to 1,3-pentadiene are compared to available C3-C5 alkenes and dienes, with external addition calculated to be kinetically favored over internal addition. However, the correlation between heats of formation and energy barriers for H atom addition to 1,2-dienes is different from that for 1,3- and 1,4-dienes. Hydrogen atom addition and abstraction rate constants are also compared for 1,3-pentadiene, with addition found to be dominant. The subsequent unimolecular reactions on the C5H9 PES are found to be highly complex with reactions taking place on a multiple-well multiple-channel PES. For clarity, the reaction mechanism and kinetics of each C5H9 radical are discussed individually in terms of the computed enthalpy of the reaction and activation, the transition-state structure/reaction class, and also in terms of the combustion species for which the reactions are of potential importance. The reactions on the C5H9 PES are divided into three reaction classes (H-shift isomerization, cycloaddition, and ß-scission reactions), and the reactivity-structure-based estimation rules for energy barriers are derived for these three reaction classes and compared to literature results for alkyl radicals.

Theoretical, Experimental, and Modeling Study of the Reaction of Hydrogen Atoms with 1- and 2-Pentene.

Alkyl radicals are prominent in combustion chemistry as they are formed by hydrocarbon decomposition or from a radical attack on hydrocarbons. Accurate determinations of the thermochemistry and kinetics of their unimolecular isomerization and decomposition reactions and related addition reactions of alkenes are therefore important in simulating the combustion chemistry of virtually all hydrocarbon fuels. In this work, a comprehensive potential energy surface (PES) for H-atom addition to and abstraction from 1- and 2-pentene, and the subsequent C-C and C-H ß-scission reactions, and H-atom transfer reactions has been considered. Thermochemical values for the species on the C5H11 PES were calculated as a function of temperature (298-2000 K), with enthalpies of formation determined using a network of isodesmic reactions. High-pressure limiting and pressure-dependent rate constants were calculated using the Rice-Ramsperger-Kassel-Marcus theory coupled with a one-dimensional master equation. As a validation of our theoretical results, hydrogen atomic resonance absorption spectrometry experiments were performed on the H-atom addition and abstraction reactions of 1- and 2-pentene. By incorporating our calculations into a detailed chemical kinetic model (AramcoMech 3.0), excellent agreement with these experiments is observed. The theoretical results are further validated via a comprehensive series of simulations of literature data. Our a priori model is found to reproduce important absolute species concentrations and product ratios reported therein.

Extensive Theoretical Study of the Thermochemical Properties of Unsaturated Hydrocarbons and Allylic and Super-Allylic Radicals: The Development and Optimization of Group Additivity Values.

In this study, the thermochemistry of C2-C7 unsaturated hydrocarbons (22 alkene and 6 diene molecules) and 16 allylic and 5 super-allylic radicals is determined using high-accuracy quantum chemistry calculations. In addition, the group additivity values (GAVs) of a total of 19 relevant groups are systematically optimized on the basis of the calculated thermochemistry of species clusters. The M06-2X method using the 6-311++G(d,p) basis set is used for the geometry optimizations, vibrational frequency calculations, and internal rotation scans for lower-frequency modes. The composite compound methods, CBS-APNO, G3, and G4, are utilized to derive the average atomization formation enthalpies. The entropy and temperature-dependent heat capacity values of all species are calculated using statistical thermodynamics in MultiWell. These results are in good agreement with literature data. A GAVs optimization is performed on the basis of a statistical analysis: a Bland-Altman plot, which is employed to visualize the agreement between the results from the quantum chemical calculations and the GA method. It is found that the 298 K entropies of the CD/C2, C/CD2/H2, C/C/CD2/H, and C/CD3/H groups disagree by more than 5 cal K-1 mol-1 compared to existing values, while the values for the ALLYLS and ALLYLT radical groups also differ by ∼2.4 and 4.1 cal K-1 mol-1, respectively. The 298 K formation enthalpies of the C/CD2/H2, C/C/CD2/H, C/CD3/H, and ALLYLT groups are modified by more than 1 kcal mol-1, compared to existing values. The updated GAVs can be used with increased confidence to estimate the thermochemical properties of combustion-relevant unsaturated hydrocarbon molecules and their radicals which are critical for the development of accurate chemical kinetic models describing the pyrolysis and oxidation of hydrocarbon and oxygenated hydrocarbon fuels.

Theoretical Kinetics Analysis for H Atom Addition to 1,3-Butadiene and Related Reactions on the C4H7 Potential Energy Surface.

The oxidation chemistry of the simplest conjugated hydrocarbon, 1,3-butadiene, can provide a first step in understanding the role of polyunsaturated hydrocarbons in combustion and, in particular, an understanding of their contribution toward soot formation. On the basis of our previous work on propene and the butene isomers (1-, 2-, and isobutene), it was found that the reaction kinetics of H-atom addition to the C═C double bond plays a significant role in fuel consumption kinetics and influences the predictions of high-temperature ignition delay times, product species concentrations, and flame speed measurements. In this study, the rate constants and thermodynamic properties for H-atom addition to 1,3-butadiene and related reactions on the C4H7 potential energy surface have been calculated using two different series of quantum chemical methods and two different kinetic codes. Excellent agreement is obtained between the two different kinetics codes. The calculated results including zero-point energies, single-point energies, rate constants, barrier heights, and thermochemistry are systematically compared among the two quantum chemical methods. 1-Methylallyl (C4H71-3) and 3-buten-1-yl (C4H71-4) radicals and C2H4 + C2H3 are found to be the most important channels and reactivity-promoting products, respectively. We calculated that terminal addition is dominant (>80%) compared to internal H-atom addition at all temperatures in the range 298-2000 K. However, this dominance decreases with increasing temperature. The calculated rate constants for the bimolecular reaction C4H6 + H → products and C2H4 + C2H3 → products are in excellent agreement with both experimental and theoretical results from the literature. For selected C4 species, the calculated thermochemical values are also in good agreement with literature data. In addition, the rate constants for H atom abstraction by H atoms have also been calculated, and it is found that abstraction from the central carbon atoms is the dominant channel (>70%) at temperatures in the range of 298-2000 K. Finally, by incorporating our calculated rate constants for both H atom addition and abstraction into our recently developed 1,3-butadiene model, we show that laminar flame speed predictions are significantly improved, emphasizing the value of this study.

Chemical Kinetics of Hydrogen Atom Abstraction from Allylic Sites by 3O2; Implications for Combustion Modeling and Simulation.

Hydrogen atom abstraction from allylic C-H bonds by molecular oxygen plays a very important role in determining the reactivity of fuel molecules having allylic hydrogen atoms. Rate constants for hydrogen atom abstraction by molecular oxygen from molecules with allylic sites have been calculated. A series of molecules with primary, secondary, tertiary, and super secondary allylic hydrogen atoms of alkene, furan, and alkylbenzene families are taken into consideration. Those molecules include propene, 2-butene, isobutene, 2-methylfuran, and toluene containing the primary allylic hydrogen atom; 1-butene, 1-pentene, 2-ethylfuran, ethylbenzene, and n-propylbenzene containing the secondary allylic hydrogen atom; 3-methyl-1-butene, 2-isopropylfuran, and isopropylbenzene containing tertiary allylic hydrogen atom; and 1-4-pentadiene containing super allylic secondary hydrogen atoms. The M06-2X/6-311++G(d,p) level of theory was used to optimize the geometries of all of the reactants, transition states, products and also the hinder rotation treatments for lower frequency modes. The G4 level of theory was used to calculate the electronic single point energies for those species to determine the 0 K barriers to reaction. Conventional transition state theory with Eckart tunnelling corrections was used to calculate the rate constants. The comparison between our calculated rate constants with the available experimental results from the literature shows good agreement for the reactions of propene and isobutene with molecular oxygen. The rate constant for toluene with O2 is about an order magnitude slower than that experimentally derived from a comprehensive model proposed by Oehlschlaeger and coauthors. The results clearly indicate the need for a more detailed investigation of the combustion kinetics of toluene oxidation and its key pyrolysis and oxidation intermediates. Despite this, our computed barriers and rate constants retain an important internal consistency. Rate constants calculated in this work have also been used in predicting the reactivity of the target fuels of 1-butene, 2-butene, isobutene, 2-methylfuran, 2,5-dimethylfuran, and toluene, and the results show that the ignition delay times for those fuels have been increased by a factor of 1.5-3. This work provides a first systematic study of one of the key initiation reaction for compounds containing allylic hydrogen atoms.

Development of a Joint Hydrogen and Syngas Combustion Mechanism Based on an Optimization Approach.

A comprehensive and hierarchical optimization of a joint hydrogen and syngas combustion mechanism has been carried out. The Kéromnès et al. (Combust Flame, 2013, 160, 995-1011) mechanism for syngas combustion was updated with our recently optimized hydrogen combustion mechanism (Varga et al., Proc Combust Inst, 2015, 35, 589-596) and optimized using a comprehensive set of direct and indirect experimental data relevant to hydrogen and syngas combustion. The collection of experimental data consisted of ignition measurements in shock tubes and rapid compression machines, burning velocity measurements, and species profiles measured using shock tubes, flow reactors, and jet-stirred reactors. The experimental conditions covered wide ranges of temperatures (800-2500 K), pressures (0.5-50 bar), equivalence ratios (ϕ = 0.3-5.0), and C/H ratios (0-3). In total, 48 Arrhenius parameters and 5 third-body collision efficiency parameters of 18 elementary reactions were optimized using these experimental data. A large number of directly measured rate coefficient values belonging to 15 of the reaction steps were also utilized. The optimization has resulted in a H2/CO combustion mechanism, which is applicable to a wide range of conditions. Moreover, new recommended rate parameters with their covariance matrix and temperature-dependent uncertainty ranges of the optimized rate coefficients are provided. The optimized mechanism was compared to 19 recent hydrogen and syngas combustion mechanisms and is shown to provide the best reproduction of the experimental data.

Modeling Nitrogen Species as Pollutants: Thermochemical Influences.

To simulate emissions of nitrogen-containing compounds in practical combustion environments, it is necessary to have accurate values for their thermochemical parameters, as well as accurate kinetic values to describe the rates of their formation and decomposition. Significant disparity is observed in the literature for the former, and we therefore present herein high-accuracy ab initio gas-phase thermochemistry for 60 nitrogenous compounds, many of which are important in the formation and consumption chemistry of NOx species. Several quantum-chemical composite methods (CBS-APNO, G3, and G4) were utilized to derive enthalpies of formation via the atomization method. Entropies and heat capacities were calculated from traditional statistical thermodynamics, with oscillators treated as anharmonic based on ro-vibrational property analyses carried out at the B3LYP/cc-pVTZ level of theory. The use of quantum chemical methods, along with the treatments of anharmonicities and hindered rotors, ensures accurate enthalpy of formation, entropy, and heat capacity values across the temperature range 298.15-3000 K. The implications of these results for atmospheric and combustion modeling are discussed.

Toward the Development of a Fundamentally Based Chemical Model for Cyclopentanone: High-Pressure-Limit Rate Constants for H Atom Abstraction and Fuel Radical Decomposition.

Theoretical aspects of the development of a chemical kinetic model for the pyrolysis and combustion of a cyclic ketone, cyclopentanone, are considered. Calculated thermodynamic and kinetic data are presented for the first time for the principal species including 2- and 3-oxo-cyclopentyl radicals, which are in reasonable agreement with the literature. These radicals can be formed via H atom abstraction reactions by H and Ö atoms and È®H, HÈ®2, and CH3 radicals, the rate constants of which have been calculated. Abstraction from the ß-hydrogen atom is the dominant process when È®H is involved, but the reverse holds true for HÈ®2 radicals. The subsequent ß-scission of the radicals formed is also determined, and it is shown that recent tunable VUV photoionization mass spectrometry experiments can be interpreted in this light. The bulk of the calculations used the composite model chemistry G4, which was benchmarked in the simplest case with a coupled cluster treatment, CCSD(T), in the complete basis set limit.

Modeling Ignition of a Heptane Isomer: Improved Thermodynamics, Reaction Pathways, Kinetics, and Rate Rule Optimizations for 2-Methylhexane.

Accurate chemical kinetic combustion models of lightly branched alkanes (e.g., 2-methylalkanes) are important to investigate the combustion behavior of real fuels. Improving the fidelity of existing kinetic models is a necessity, as new experiments and advanced theories show inaccuracies in certain portions of the models. This study focuses on updating thermodynamic data and the kinetic reaction mechanism for a gasoline surrogate component, 2-methylhexane, based on recently published thermodynamic group values and rate rules derived from quantum calculations and experiments. Alternative pathways for the isomerization of peroxy-alkylhydroperoxide (OOQOOH) radicals are also investigated. The effects of these updates are compared against new high-pressure shock tube and rapid compression machine ignition delay measurements. It is shown that rate constant modifications are required to improve agreement between kinetic modeling simulations and experimental data. We further demonstrate the ability to optimize the kinetic model using both manual and automated techniques for rate parameter tunings to improve agreement with the measured ignition delay time data. Finally, additional low temperature chain branching reaction pathways are shown to improve the model's performance. The present approach to model development provides better performance across extended operating conditions while also strengthening the fundamental basis of the model.

Experimental and Kinetic Modeling Study of 2-Methyl-2-Butene: Allylic Hydrocarbon Kinetics.

Two experimental studies have been carried out on the oxidation of 2-methyl-2-butene, one measuring ignition delay times behind reflected shock waves in a stainless steel shock tube, and the other measuring fuel, intermediate, and product species mole fractions in a jet-stirred reactor (JSR). The shock tube ignition experiments were carried out at three different pressures, approximately 1.7, 11.2, and 31 atm, and at each pressure, fuel-lean (ϕ = 0.5), stoichiometric (ϕ = 1.0), and fuel-rich (ϕ = 2.0) mixtures were examined, with each fuel/oxygen mixture diluted in 99% Ar, for initial postshock temperatures between 1330 and 1730 K. The JSR experiments were performed at nearly atmospheric pressure (800 Torr), with stoichiometric fuel/oxygen mixtures with 0.01 mole fraction of 2M2B fuel, a residence time in the reactor of 1.5 s, and mole fractions of 36 different chemical species were measured over a temperature range from 600 to 1150 K. These JSR experiments represent the first such study reporting detailed species measurements for an unsaturated, branched hydrocarbon fuel larger than iso-butene. A detailed chemical kinetic reaction mechanism was developed to study the important reaction pathways in these experiments, with particular attention on the role played by allylic C-H bonds and allylic pentenyl radicals. The results show that, at high temperatures, this olefinic fuel reacts rapidly, similar to related alkane fuels, but the pronounced thermal stability of the allylic pentenyl species inhibits low temperature reactivity, so 2M2B does not produce "cool flames" or negative temperature coefficient behavior. The connections between olefin hydrocarbon fuels, resulting allylic fuel radicals, the resulting lack of low-temperature reactivity, and the gasoline engine concept of octane sensitivity are discussed.

Revisiting the Kinetics and Thermodynamics of the Low-Temperature Oxidation Pathways of Alkanes: A Case Study of the Three Pentane Isomers.

This paper describes our developing understanding of low-temperature oxidation kinetics. We have investigated the ignition of the three pentane isomers in a rapid compression machine over a wide range of temperatures and pressures, including conditions of negative temperature coefficient behavior. The pentane isomers are small alkanes, yet have structures that are complex enough to allow for the application of their kinetic and thermochemical rules to larger molecules. Updates to the thermochemistry of the species important in the low-temperature oxidation of hydrocarbons have been made based on a thorough literature review. An evaluation of recent quantum-chemically derived rate coefficients from the literature pertinent to important low-temperature oxidation reaction classes has been performed, and new rate rules are recommended for these classes. Several reaction classes have also been included to determine their importance with regard to simulation results, and we have found that they should be included when developing future chemical kinetic mechanisms. A comparison of the model simulations with pressure-time histories from experiments in a rapid compression machine shows very good agreement for both ignition delay time and pressure rise for both the first- and second-stage ignition events. We show that revisions to both the thermochemistry and the kinetics are required in order to replicate experiments well. A broader validation of the models with ignition delay times from shock tubes and a rapid compression machine is presented in an accompanying paper. The results of this study enhance our understanding of the combustion of straight- and branched-chained alkanes.

Theoretical chemical kinetic study of the H-atom abstraction reactions from aldehydes and acids by H atoms and ȮH, HȮ2, and CH3 radicals.

We have performed a systematic, theoretical chemical kinetic investigation of H atom abstraction by H atoms and ȮH, HȮ2, and CH3 radicals from aldehydes (methanal, ethanal, propanal, and isobutanal) and acids (methanoic acid, ethanoic acid, propanoic acid, and isobutanoic acid). The geometry optimizations and frequencies of all of the species in the reaction mechanisms of the title reactions were calculated using the MP2 method and the 6-311G(d,p) basis set. The one-dimensional hindered rotor treatment for reactants and transition states and the intrinsic reaction coordinate calculations were also determined at the MP2/6-311G(d,p) level of theory. For the reactions of methanal and methanoic acid with H atoms and ȮH, HȮ2, and CH3 radicals, the calculated relative electronic energies were obtained with the CCSD(T)/cc-pVXZ (where X = D, T, and Q) method and were extrapolated to the complete basis set limit. The electronic energies obtained with the CCSD(T)/cc-pVTZ method were benchmarked against the CCSD(T)/CBS energies and were found to be within 1 kcal mol(-1) of one another. Thus, the energies calculated using the less expensive CCSD(T)/cc-pVTZ method were used in all of the reaction mechanisms and in calculating our high-pressure limit rate constants for the title reactions. Rate constants were calculated using conventional transition state theory with an asymmetric Eckart tunneling correction, as implemented in Variflex. Herein, we report the individual and average rate constants, on a per H atom basis, and total rate constants in the temperature range 500-2000 K. We have compared some of our rate constant results to available experimental and theoretical data, and our results are generally in good agreement.